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Abstract 

A review on the current efforts to approach and to surpass the fundamental limit in 
the sensitivity of the Weber type gravitational wave antennae is reported. Applications of 
quantum non- demolition techniques to the concrete example of an antenna resonant with 
the transducer are discussed in detail. Analogies and differences from the framework of the 
squeezed states in quantum optics are finally discussed. 


1 Introduction 

The importance of detecting gravitational waves, as frequently pointed out, consists not only 
in verifying one of the most direct and astonishing predictions of the simplest metric theory of 
gravitation, i.e. General Relativity, but also in the possibility to open new windows on phenomena 
in the Universe in which only violent releases of gravitational energy occur [1]. Gravitational waves 
have not yet been directly measured because of the extreme smallness of the power emitted even 
by astronomical systems. The hypothetical sources which are strong candidates for emitting 
gravitational waves, according to our understanding of them actually only due to informations 
collected via the electromagnetic astronomy, are divided into two classes based upon the time 
evolution. Impulsive sources can be catastrofic events such as supernovae explosions and collapsing 
binary systems. The frequency spectrum of gravitational waves of this kind is flat up to 10 3 Hz, 
these impulsive phenomena having a characteristic duration of the order of milliseconds. One 
expects a perturbation of the metric tensor h ~ 10~ 21 — 10“ 18 for events in our Galaxy and 
h % 10- 23 — 10- 21 for events in the Virgo Cluster. Periodic sources can be pulsars if they deviate 
substantially from axial symmetry. The expected frequencies range is in this case between 10 -2 and 
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10 2 Hz, while h % 10“ 27 — 10 -25 . The efforts to detect gravitational waves have been concentrated 
from the very beginning on the impulsive events because of the larger expected perturbation to 
the metric tensor. It turns out that the modulation of the space-time induced by a gravitational 
wave on an extended body can also be seen as a production of a force field in it. Detecting 
the gravitational wave is therefore translated into the problem of detecting this small force of 
geometrical nature and the displacements produced by it in a test mass. The displacement induced 
in a body of reasonable sizes, %lm, has therefore an amplitude of the order of 10“ 21 if the event 
is due to the a supernovae event in the Virgo Cluster. The accuracy required to measure such a 
small displacement is so high that the quantum nature of the detector has to be taken into account 
because the De Broglie wavelenght of a macroscopic test mass is of the same order of magnitude of 
the expected signal due to the gravitational waves. Here we report on the status of the art of the 
measurement techniques developed to allow monitoring of a class of gravitational wave detectors 
in a quantum regime. After a brief introduction for schematizing the detectors of gravitational 
waves and the sensitivity limit due to the fundamental noise in part 2, we introduce, in part 3, the 
quantum non-demolition measurement schemes for overcoming these limitations. The applications 
of stroboscopic and continuous quantum non-demolition schemes for a gravitational bar antenna 
resonant with the transducer are described respectively in part 4 and part 5. Conclusions deal also 
with the analogies and the differences from the quantum optics framework and the importance of 
this topic for understanding quantum mechanics applied to single macroscopic degrees of freedom 
repeatedly monitored. 

2 Weber gravitational antennae: fundamental sensitivity 
limits 

The gravitational wave detectors devised so far are based upon monitoring of the distance between 
two masses localized at different points. The equivalence principle requires a non-local, extended, 
structure of a gravitational wave detector because it is possible to nullify locally the effects of a 
gravitational field by means of a suitable choice of the reference frame. 

Let us consider two masses in free fall: what is then measured is their variable distance which is 
supposed to be much smaller than the gravitational wavelength. The effect of a gravitational wave 
coming along z axis with proper polarization is to increase of hj 2 the distance along y axis and 
to decrease by h/2 the distance along x axis. A classification of the gravitational wave detectors 
divides these into non resonant and resonant detectors if the two masses are respectively free or 
elastically coupled. 

In non resonant detectors the distance between the two masses is measured by means of 
interferometric devices. The arms of the interferometer proposed so far are of the order of Km and 
use of multiple reflections allows to increase the physical paths by several orders of magnitude. 
In this contribution we will not be concerned with this kind of detectors but we shall instead 
consider the resonant detectors (Weber type gravitational wave antennae), the quantum limit in a 
interferometric antenna being enforced by the shot noise and the momentum fluctuations imparted 
by the photon flux to the central mirrow of the interferometer [2]. 

Resonant antennae are tipically cylindrical bars of materials having low internal dissipations. 
The materials used are silicon, sapphire, niobium or a particular aluminum alloy (A1 5056) and 
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the mass of the antennae is of few tons. 

One can show that the motion of the ends of a cylindrical bar of mass M and length L 
oscillating in its fundamental longitudinal mode is equivalent to that of a harmonic oscillator of 
mass M/2 and equivalent length 4 L/tt 2 . If x is the displacement from equilibrium position the 
equation of motion of the Weber oscillator is 


i + — +u >Zx = -^-Lh(t) (1) 

To 7T^ 

where r 0 is the damping time, u> 0 is the proper frequency and &(t) is the amplitude of the incoming 
gravitational wave. The forcing term due to the gravitational field is proportional to the distance 
between the two masses. From this formula one can calculate the cross section for the transfer 
of energy from the wave to the antenna and one finds that this is proportional to the mass of 
the antenna and to L 2 . The proper frequency is chosen to be tuned with the frequency of the 
expected wave (10 3 Hz) and the corresponding wavelenght is very large compared to the size of the 
antenna. To amplify the extremely small oscillations coupling of the bar with another oscillator 
of very small mass is used [3], [4]. In this case a system of two coupled harmonic oscillators is 
obtained in which the energy is continuosly transferred back and forth from M to m via beating. 
If the dissipations in the two oscillators are made negligible the amplitude of the oscillations in 
the second resonator is increased by a factor 1 /y/Ji with respect to the first resonator, where 
p = m/M, provided that the frequencies of the two uncoupled oscillators are made coincident. 
The motion of the transducer is transformed into an electric signal by means of a variable capacity 
and an amplifier schematizable as an ideal amplifier of gain A and two noise sources generators 
with current and voltage spectral densities respectively SI n and SV n . The sources of noise are 
the thermal noise, i.e. Brownian motion of antenna, which gives a contribution KT to the energy 
of the oscillator, being K the Boltzmann’s constant and T the thermodynamical temperature 
of the antenna and the amplifier noise, which is expressed by means of the parameter T n — 
(SV n Sin) 1 / 2 / Kb) called noise temperature of the amplifier. This last noise has two effects: it 
contributes directly as an additive noise source at the output and it acts on the transducer leading 
to an increase of the temperature. In other words every transducer is at the same time an actuator 
and the amplifier noise gives rise to a back-action force acting on the mechanical oscillator. 

If we define a noise temperature T e ff as the temperature which corresponds to the minimum 
detectable energy E e ff = KT e ff transferred to the bar by an impulsive signal with an output 
signal/noise ratio equal to 1, we find, using a Wiener algorithm in the data analysis [5] 



( 2 ) 


where Q = u> 0 t 0 is the quality factor of the mechanical system, /3 is substantially the fraction of 
energy transferred to the electromagnetic circuit by the bar throug the capacitive coupling and 
A 0 the impedance matching factor defined as 


= SV n 1 
SI n Zq' 


( 3 ) 


For the antenna of the Rome group continously operating since one year at CERN one has a 
thermodynamical temperature of ^ 4.2K; the other parameters are Q ~ 5 ■ 10 6 and an amplifier 
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noise temperature T n ~ 10“ 7 K [6]. It has been possible to achieve this last result making use of 
a SQUID amplifier. So one gets for T e // a value of ^ lOfiK, which is not far from the quantum 
limit temperature 

T qL = Y~ 1(r * K (4) 

One expects that the force with which a gravitational wave acts on the antenna is by many 
orders of magnitude below the thermal noise even at thermodynamical temperatures as low as 
lOmK which is the temperature at which the third generation antennae will operate. Howewer, 
due to the particular features of the data analysis based on the variation of energy in the oscillator 
in the time, the quantum regime is reached earlier than as expected by (4). By writing the amount 
of energy which is exchanged during the measurement time At between the harmonic oscillator 
and the thermal reservoir and the quantized energy introduced by the measuring apparatus is 
easy to show that the quantum regime is obtained when the following condition is satisfied 


KTAt 


<h 


( 5 ) 


This can be also shown by reasoning in terms of displacements instead of energy. The variation 
of the length of the bar due to a gravitational wave with amplitude h is, according to (1) 


M 

l 



( 6 ) 


Because typical values for h are h = 10“ 21 (which corresponds to a supernova explosion in 
the center of the Galaxy) taking I = lm, one gets from (5) a variation of the length of the bar 
AL ~ 10“ 19 cm which coincides with the standard quantum limit (i.e. the root square mean of 
the position of a harmonic oscillator in his fundamental mode) 


< 7) 

It follows therefore that if we do not overcome this limit no information can be obtained on the 
evolution of the harmonic oscillator. 

In these conditions one can find a method to measure the position of the quantum oscillator 
and to see if an external force has acted on it. However in doing this one must take into account 
that the position operator x(t) does not commute with itself at different times. Indeed with a 
measurement of x(t) at time t one put the oscillator into an eigenstate of x(t); if one repeats this 
measurement at the instant t + r one puts the oscillator into another eigenstate. It turns out 
that it is not possible to know if the change in x(t) is caused by a very weak classical external 
force because of the quantum demolition of the state. What is needed is therefore a measurement 
which does not prevent the execution of the next measurements of the same observable avoiding 
the demolition of the projection of the state on that observable. This is possible in non-relativistic 
quantum mechanics as we will discuss in the following considerations, because this theory make 
limitations only on a simultaneous, perfect knowledge of two canonical observables. 
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3 Quantum non demolition measurements 

The introduction of the quantum non-demolition measurements (QND) dates back to an article 
by Landau and Peierls [7] in 1931. However only recently, after understanding the role of quantum 
mechanics in the fundamental limits to the amplifier sensitivity [8], [9] and under the request to 
surpass the quantum limit in detectors of small displacements [10], [11], the problem has been 
studied in detail [12], [13]. The idea of a QND strategy is to perform a series of measurements 
of one observable of a single object in such a way that the act of the measurement itself does 
not affect the predictability of the result of the next measurements of the same observable. In 
order to do this the observable, the instants of time in which it is observed and the interaction 
H amil tonian should be all carefully chosen for a given dynamical system. For instance, a first high 
precision measurement of the position of a free particle implies a large dispersion in the possible 
values of measurements of momentum. If a second measurement of position is made, due to the 
Heisenberg evolution, the result will have a large dispersion too. Instead, if a measurement of 
momentum in a free particle is made at a given instant of time, a second possible measurement 
will give the same result due to the constant value of the momentum between the two consecutive 
measurement, provided that the interaction dtie to the first measurement has not demolished the 
state. This simple example shows the route to define quantum non-demolition measurements. 
Only particular observables which satisfy a commutation relation at different times U and tj are 
allowed to be monitored in a QND way, i.e. if 

[*(*),*(«,-)] = 0 . ( 8 ) 

Moreover, we must also take into account the perturbation on £(£) induced by the measuring 
apparatus which is coupled to the observed system by means of the operator Hamiltonian Hi . 
To avoid changes in the expected value of the observable during the measurement the following 
condition must be satisfied: 


[*(*),£] = 0. (9) 

This condition assures that the interaction Hamiltonian is simultaneously diagonalizable with the 
measured observable, no changes are induced in the measured observable during the measurement 
time in which only the interaction Hamiltonian will be responsible for the time evolution. A 
sequence of measurements performed under conditions (8) and (9) will give always the same result. 
This is a definition of a QND measurement. If the instants of time in which it is satisfied (8) are 
discrete the QND scheme is named stroboscopic or, in a realistic configuration with a duration of 
the measurement small with respect to the characteristic timescale of the motion of the observed 
system, quasi-stroboscopic [14], [15], [16]. Otherwise, having a continuous set of instants of time, 
the QND scheme is named continuous. 

In the case of a single oscillator one introduces the two components of the complex amplitude 

; do 

X 2 = Im[(x + i-^ i )e iut } 

such that x{t) = X 2 cos ut + X 2 sin u>t. Their properties are 
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(a) 

dX x dX 2 _ 
dt dt 

0 [A"i(t), X x (t + r)| = [X 2 (0> X 2 (t + r)] = 0 

(ii) 


(b) 

[X l (0,Xi(<)] = — • 

mw 

(12) 


By using (a) and (6) we get 


[i(t),i(t + r)] = ~[Xi,X 2 ] {cos u>t sin u)(t + r) — sin wt cos u>(t + r)} = — — sin u/r. (13) 

mu> 

This means that to do a QND measurement of the operator i(t) in a single harmonic oscillator 
one needs the Hamiltonian (here q is the variable of the measuring apparatus which couples with 
the oscillator) 


Hi i = E 0 6(t )xq (14) 

U) 

such that the interaction between the system and the measuring apparatus is turn on only when 
£(*) commutes with itself, that is why this kinds of measurements are called stroboscopic Q.N.D.. 
For a component of the complex amplitude, X \ , a QND interaction Hamiltonian should be 

[ 12 ] 


Hi — EqX\<i 


that is approximately obtained by using the interaction Hamiltonian 


(15) 


Hi — 2Eo cos u> m t x q 


(16) 


provided a low-pass filter at w c << u> m is used. For practical reasons a different pumping is used, 
namely a up-conversion around an electrical frequency u> e such that the interaction Hamiltonian 
is now 


A 2£ 

Hi = E 0 cos u > € t cos uj m t x q = — [cos(u; e -f u; m )t + cos(w e — w m )£] x q 

z 


(17) 


which allows an approximate measurement of X\ if a filtering around o> e is performed with a 
selectivity such that the terms oscillating at ± 2u/ m are made negligible. It has been pointed 
out that the continous approximate QND measurement scheme of one component of the complex 
amplitude is obtained as a first order approximation of the corresponding stroboscopic scheme 
[17]. If we start from the interaction Hamiltonian of a stroboscopic measurement of X\ expressed 
in terms of the physical observable x 


Hi = E 0 cos u> e t £ 6(t - — )X x q = E 0 cos u; e t £(- l) n ^(i - — ) x q 

n W 1 n 

we will see that, by Fourier expanding the Dirac-distribution, it is obtained 

Hi = E 0 cos cos(2n + l)u»i tx q 


(18) 


(19) 
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that, at the first order, is 


Hi — E 0 cos t cos cjjtx q (20) 

i.e. the usual approximate scheme for monitoring of X \ . Thus knowing a QND stroboscopic 
strategy it is simple to write the corresponding QND approximate continuous strategy. This 
property will be particularly useful in the following considerations, where the more complicated 
but realistic case of two coupled harmonic oscillators will be treated. 

It has been pointed out that also in the classical regime, i.e. when the amplifier is not quantum 
limited, the QND measurement schemes provide a better sensitivity because one phase of the signal 
is shielded by the back-action force of the amplifier. A quantitative model in the classical limit 
has been developed in [18]: it turns out that by writing the noise temperature as 

T b = — T n - (21) 

U> e T 


for a standard ’amplitude and phase’ monitoring is r < 1, and for a QND/BAE scheme r may 
be greater than unity. This is due to the squeezing of the electrical noise into one mechanical 
phase. A generalized uncertainty relation for the two classical conjugate observables due to the 
back-action of the amplifier noise is introduced as 


AX x AX 2 


K B T n 
2 mu? m uj e 


( 22 ) 


which may be obtained through a replacement on the right hand side in the standard quantum 
uncertainty relationship 

AX^AXj ~ — — (23) 

2mw ra v ' 

of h with K B T n j(jj 2 . If a squeezing factor p such that AX\ = pAX 2 is introduced (p — > 0 means 
a noise-free measurement of X x ) the minimum burst noise temperature can be written as 


Tb 


rruJL AX} 1 _ w\ 

~ -T n —p 

2 4 u>2 


(24) 


showing that the r figure of merit has a dynamical interpretation in terms of a squeezing factor. 
Recently, an interpretation of the back-action evasion strategies in which they are seen as an 
alternative to the usual impedance matching for maximizing the signal to noise ratio has been 
discussed [19]. 

The description of the QND measurement suggests how to measure small forces below the 
standard quantum limit. By means of a simple integration of the Heisenberg equation in presence 
of an external force F(t), one gets for the QND operator X\ 


Xiit) = X 1 (t 0 ) - 1 t — sin wt'dt'. (25) 

Jt 0 TTlijJ 

A sequence of measurement of X x will then give as a result a sequence of eigenstates linked to the 
value of the external force 
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{(«,r) = tft 0 )- [' ™ An wt 9 ,*'. ( 26 ) 

v Jto mu 

By means of successive measurements it is possible to study the form of F(t) simply inverting ( 26 ) 


771a; d 




( 27 ) 


The singularities for t = n n/u corresponds to a null information on the force acting on the 
harmonic oscillator on some instants of time. This can be compensated by using a second oscillator 
(i.e. a second antenna) with complex amplitude Y\ + iYj which has eigenvalues 


C(*,r) = C(M- — «ia 

J tg TTUjJ 

here obviously the singularities are in t n = ( 2 n + l)w/ 2 u). 


( 28 ) 


4 QND quasi-stroboscopic scheme for coupled harmonic 
oscillators 

The current generation of gravitational wave antenna of the Weber type operates by means of 
an antenna coupled to a small mechanical resonator. In such a way the energy deposited in the 
antenna by a gravitational wave burst is transferred to the transducer. In the case of an ideal 
transfer of energy, i.e. with both a perfect tuning of the two uncoupled frequencies and negligible 
dissipations during the beating period, the amplitude of the oscillations in the transducer is larger 
than that in the antenna by a factor equal to the square root of the ratio of the equivalent masses 
of the two resonators. All the detectors operating in coincidence as described in [ 6 ] were equipped 
with a resonant transducer and the same is also planned for the third generation of gravitational 
wave antennas cooled at 50 mK now under development. It is therefore important to generalize 
the previous considerations on the QND schemes to this situation, as already outlined in [ 20 ]. As 
we have seen, it is possible to schematize the gravitational cryogenic antenna and the resonant 
transducer with two coupled harmonic oscillator having masses respectively m x and m v (with 
H — IJbt 1 ). The two coupled mechanicals oscillators are described by the Lagrangian 


L = ^rn m x 2 + ^rn v y 2 - x 3 - ^m y u y 3 {y - i) 1 = 17) T 


p9) 


where the normalized coordinates £ = y/m m x and 77 = y/rn^y have been introduced, together with 
the matrices T and V 
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As we have already cited to obtain the maximum coupling the two oscillators should have the 
same frequency w* = ut y — u;, i.e. they should be tuned. In this case one finds the solutions 


u/ ± 


="o’(i+f±\H i +j)) 


(32) 


which we can write more easily introducing a± = ^ ^1 + obtaining — u/q(1 + ot±). 

The normal coordinates corresponding to the eigenfrequendes u)± are linked to the physical 
coordinates by means of an orthogonal matrix 


(lr)- 


I- 7-1 ?= ' 

ya_ + 2 vW a - ^ 

1 a 

V \/ a + + 2 v / M\/ a+ ” ^ f 




(33) 


Let us introduce the complex amplitudes of the normal modes 

r ± _ A . PS 


which satisfy the relations 


as well as 


Xi = E± cosa >±t — sinu;±< 

^ p_ 

.Yf = E± sinu>±t + cos w±t 


[Xt,X}) = 2- [X;,X; 1 = £ 

UJ+ u>_ 


(34) 


[Y+ 2 (t),jr+ 2 (« + r)] = [X^(t) t Xr A (t + r)] = 0. 

We can also rewrite the Hamiltonian H of the system as 

h = ^{(x+r + (x+Y ] + y [(xn 2 + (x 2 -) 3 }- 


(35) 


(36) 


(37) 


The commutator [jK*X“ x(t) y y(t + r) — x(t + r)] is calculated by writing y and x in terms of the 
complex amplitudes °f the normal modes which are integral of the motion and by using 

the same computation procedure which led us to formula (13). Using (35), (36) we obtain, finally, 
the expression 


[y(0 - *{t),y{t + t)~ x(t + r)] = 


ih 


+ 4 




to* 


u>Z 


sm u ;+r H sin u ;_r 


This quantity becomes, in the limit \i — ► 0 

[y(t) - x(t),y{t + r) - x(t + r)] = 


ih <jj\ + u) 2 _ 

2 m y u) u) 2 




sin lot cosubt 


(38) 


(39) 
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where u> = + 4 — ► and a;# = — 1 g — = f \/^* The (38) and (39) show that 

the commutator of the operator y — £ with itself at different times is time dependent and it 
has a characteristic beating behaviour. We have seen that in a quasi-stroboscopic scheme for 
a single harmonic oscillator the commutator is zero each half a period of the motion and the 
stroboscopicity is defined whenever measurements with a duration small compared to the period 
of the motion are performed. This implies a measurement time, a duty cycle, very small and a 
consequent small value of the effective electromechanical quality factor. In the case of a double 
harmonic oscillator this drawback is less pronounced because the commutativity is assured every 
half of a beating period for a time of the order of a period of oscillation. Thus quasi-stroboscopic 
QND schemes already proposed as a generalization of the conventional BAE scheme based upon 
a continuous monitoring [17] and already tested on a single oscillator system [21] can be adapted 
to this situation. In the case of a single harmonic oscillator the duration of the measurement must 
be small compared to the period of the harmonic oscillator T, in the case of two coupled harmonic 
oscillators this duration is of the order of some periods of the uncoupled oscillator, although the 
interaction must be turned on every quarter of a beating period. The interaction Hamiltonian for 
a two coupled harmonic oscillator system is therefore 


% = # £«' - 


nTs AT . uTb . AT X1/ . 

+ ~2~) + H-t + -$- + -2")l(y - *)« 


(40) 


where T B is the beat period and AT is of the order of the period of a single harmonic oscillator. 
Practical values are Tb — 40ms and AT ~ 2ms. To calculate the error in a quasi stroboscopic 
measurement of the operator y — x performed for instance in the interval — ^ > 3 ^ + ^7 we 
identify the conjugate observable of ij — z as the quantity ( p y — p*)/2. This last can be expressed in 
terms of the components of the amplitudes of the normal modes and the commutator at different 
times of the two conjugate observables is obtained as 

m + r) - «(< + r),i(w<) - p.m = cos ^ + 7777725 cos w - r) - (41) 

When r = 0 the commutator relationship (41) is written as 


[»(<) - -*.«)] = « ( 42 ) 

which is exactly the quantity [*(<), + [i/(^), 

By expressing u)+ and u>_ in terms of the frequencies u> and w B and substituting in a ± their 
expressions in terms of /z we get finally 


[y(f + r) - x(t + r) y -(p v (t) - p x (t))] = ih( cos ur cos t c B r sin ur sin a/ £ r) (43) 

2 \W + 4 ) 

If the measurement is performed in the interval [ 3 — — + ^], we can approximate 

C 03 u> b t ~ 1 and sin u> B r ~ u ) B r — | and a measurement of infinitesimal duration V performed 
in such interval and with a precision A[y(t) — x(t)] allows 1° evaluate the error introduced in the 
measurement process on the uncertainty product as 


218 


A[y(< + t ') - x(t + <')] • ^A[p„(f) - p*(t)] ~ ^|cos u it' 1 - + - sin ut'(u> B t' - ^)| (44) 

z z v/Mm + 4) 1 

from which, under the approximation for the trigonometric functions, we obtain 

-A\p v (t) - p x (t)} « 2A[y(0 -*(<)]• (45) 

The error due to a measurement of duration t 7 on the operator y — x is calculated starting from 
A[y(t) — i(t)] because 

A[y(< + f') - x(t + «')] Si A[y(f) - x(t)]| cos u>t' 7 === sin u)t'(u> B t' - -)| (46) 

VW + 4) 2 

If the notation now is changed defining A ( = A[p(f) — x(t)] we have 


~~ [cos ut' - — . = + it sin u)t'\w B t' — ^|) 

1 yW + 4) 2 


and in the limit of t' — ► 0 we get 


dA t 


1 + fl TTU) 


dt v^ + 4 ) 2 


A t 


(47) 


(48) 


from which, by integrating, we obtain the error on a measurement performed around t = as 

Zujb 


A[p(f + r) - x(t + r)] Si A[p(t) - x(t)] exp[— ^ u r]. 

2y /*(/* + 4) 

Kor instance, for a choice t = and t — we obtain 


(49) 


at.-./ w , 2ir \ -/ * , 2,r M A r »/ * 2t n IT 2fl\, r 27T 2 (/l + l), 

A(y(- 1- — ) - x(- 1 ) Si A y ) - x( ) exp — = >ll_ — 1 

2<m b u> 2u> b u; 2lj b u> 2u b u + 4) 


(50) 


A drawback of these measurement scheme appears when fi is very small and the frequency of the 
measurement is consequently very small too. To overcome this problem a multimode configuration 
can be used. In this case the commutator at different times more frequently approaches zero when 
compared to a two-mode configuration of the same final mass ratio. A more detailed description 
of this point can be found in [22]. 


5 QND continuous schemes for coupled harmonic oscil- 
lators 

Also QND continuous schemes can be used for coupled harmonic oscillator. A first example is 
given by a monitoring of the complex amplitude of the physical modes x and y [23]. Introducing 
the complex amplitudes such that 
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(51) 


f y = Re[{Y x + tF 2 )e--*‘] 

l Pyl m y u) y = Img[{9 i + ty 2 )e---‘] 

we can rewrite the Hamiltonian in terms of Y x and Yj and, by writing the Heisenberg equations 
for the time evolution of Y x , we obtain 


dY x 

dt 


= —u) y x sin u/„t. 


(52) 


The complex amplitude is not a constant of the motion. However it is easily proved that it is a 
QND observable. A relationship valid for an infinitesimal time r is derived for the time evolution 


Yi(i + r ) = Y x (t) - u v x sin Wytr 
and this implies the commutation rule for Y x at different times 


(53) 


[Yi(< + r),K(0] = [KW- w » isin = 0 ( 54 ) 

because of the commutativity between Y x and x. Thus Yi (or f a , for which similar relationships 
hold) is a QND observable, although it is not conserved during the motion. From (52) the 
coordinate x is inferred as 


x(t) = - 


u v 


_1 dYi 

sin a )yt dt 


(55) 


apart from the singularities already discussed appearing when sin a > v t = 0. When a classical force 
F(t) acts on the system the Hamiltonian operator is modified and the added term is 


H, = -(* + ») F(t) 

obtaining, in this case, the following expression for the time evolution of Y x 


(56) 


dY 

dt 


—U)yi sin UJyt — 


TTlyWy 


(57) 


However the effect of the external force to be detected, in our case of geometrical nature, on 
the transducer is negligible compared to the effect on the antenna, due to the smaller size of the 
transducer. Thus H f ~ -xF’(t) and the second term in (57) can be omitted. In this reasonable 
approximation, i.e. F(t) acting only on the antenna, Y x U also QNDF, i.e. QND also in presence of 
an external force. To obtain a continuous monitoring of Y x we need a QND interaction Hamiltonian 
of the type 


Hi = E 0 cos t cos ujy t (y — x)q (58) 

that is a coherent superposition of pumpings at frequencies u ± Analogous considerations 
can be made for the monitoring of the real or the imaginary part of the complex amplitude of 
one normal mode expressed in terms of the physical modes through (33). The advantage in this 
case is that the quantity X[ + is a constant of the motion and its monitoring is the standard one 
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slrciidy discussed for & single hsnnonic oscill&tor. This is obtained by rnenns of tbe interaction 
Hamiltonian 


Hi = E a cos u> e t cos ui+ t(y — x)q (59) 

and the analogous for monitoring a component of the complex amplitude X i by substituting t*>+ 
with u>_. One drawback of monitoring one component of the complex amplitude of the normal 
modes is that the information on the other mode is lost, and it is crucial to have informations on 
both the modes to take full advantage of the resonant schemes. 

An alternative scheme suggested by the time dependence of the commutator consists in a 
monitoring corresponding to the following Hamiltonian: 


Hi = E o cos u> e t cos w t cos u>b t (y — i)q. (60) 

This coupling allows to infer informations on both the modes because, upon filtering around c j € in 
such a way to neglect terms oscillating at u> e ± 2 u>b, ± 2a>, ± 2(u> ± u > b), it can be rewritten 

as 

Hi = cos w. t{0 + X+ + 0-X;)q (61) 

4 

where j3± are coefficients related to the coefficients of the matrix (33) and are expressed as 


0 ± = [™*(2 + ^ T + J))] 1/J ( 




(±±31 

A* 


, 3 / 2 ' 


(62) 


which, in the limit of p -* 0, goes to /3 ± = In this limit the interaction Hamiltonian 

assumes a simple form 

Hi = ..fi- cos w. t(X+ - X{)q (63) 

Y O&TThxflr 

which contains informations on both the normal modes and in such a way that QND measure- 
ments can be performed on both the modes. In all the three cases here discussed the selectivity 
requirements on the electrical circuit are more stringent than in the case of a single harmonic os- 
cillator, because now the electrical oscillator must have a quality factor Q e >> u> e /u >b order to 

avoid detection of sidebands contributions. The interaction Hamiltonian (60) can also be written 
as 


Hi 


Eo 

2 


cos u> € t(cos t + cos cu- t)(y — i)q. 


(64) 


With the analogy to the multipump scheme discussed for a single oscillator we can imagine a 
interaction Hamiltonian of which (64) is only the first order approximation 


rv +00 +00 

Hi ~ COS U) e cos (2n + l)u;+ t + cos (2m + l)u/_ t](y — x)q 

2 n=0 


m— 0 


which corresponds, in the limit of a stroboscopic pumping of the kind 


+oo 


Hi = E 0 (E(-1 r*(t - — ) + E - ^)](y - *)9 

n=0 W U 


+oo 


717T, 


m=0 


(65) 


( 66 ) 
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It is interesting to observe that after a time equal to T B j 2 both the trains of Dirac distributions 
will coincide, i.e. T B f 2 = mr/u;+ — Trnr/u>- where n = m + 2 (the fact that n and m have the 
same parity assures the same sign of the corresponding Dirac pulses at those times). So each half 
a period the two trains are summed and the quasi-stroboscopic scheme discussed in the previous 
section can be considered as the first order approximation of the stroboscopic scheme resulting 
from (66). This completes the connection between the multipump continuous schemes and the 
quasi-stroboscopic scheme introduced in the previous section. 


6 Conclusions 

We have shown the scenario under which quantum non-demolition measurement schemes should 
be demanded for detecting gravitational waves in the generation of resonant gravitational wave 
antennae currently under development, particularly ultra-low temperature resonant bar antennae 
such as the Rome, Legnaro and Stanford ones which will work at a thermodynamical temperature 
of ~ 50 mK. Both QND stroboscopic and continuous schemes have been discussed as well as their 
link and practical schemes to implement them. However the interest of quantum non-demolition 
measurement schemes goes beyond the only detectability of the gravitational radiation, involving 
also the quantum measurement theory and the predictions of it for repeated measurements on 
a single macroscopic oscillator. Feasibility of the generation of macroscopically distinguishable 
states using a QND scheme has been recently discussed in quantum optics [24], [25]. It has been 
pointed out that the generation of Schroedinger cats using micromechanical oscillators with quan- 
tum limited sensitivity is also feasible [26]. Unlike the optical case, in which the QND measurement 
is obtained with a frequency mixing due to non-linear susceptivity, the QND measurement for the 
mechanical case is obtained using an electric field which can be large as one wants. Dissipations 
in a mechanical oscillator also are quite low compared to electrical or optical oscillators. More- 
over, analogies to the production and the detection of squeezed states in optics [27] have been 
shown. We want to point out a fundamental difference between the two topics: in the case of the 
optical squeezed states we deal with a quantized field in which its quantum nature is responsible 
for the limitation to the sensitivity, in the case of quantum non-demolition measurements on a 
harmonic oscillator the eventual force field which has to be monitored is considered classical and 
the fundamental limitations comes from the process of the measurement and the interaction of 
the meter with the external environment. What is squeezed in a QND measure is the back-action 
noise generated by the amplifier and the squeezing is made in a phase orthogonal to the one which 
is detected [21], Despite of this conceptual difference the formalisms to deal with QND strategies 
are similar to the one used to deal with squeezed states. This analogy is so narrow that also multi- 
pump [28], [29] and quasi-stroboscopic [30], [31] schemes have been indipendently and successfully 
implemented for squeezing the light. Further thoughts on the analogies and the differences be- 
tween quantum non-demolition measurements on a harmonic oscillator and the squeezing of the 
quantum noise can give rise to a better understanding on the same interpretation of Quantum 
Electrodynamics and the operative origin of the vacuum fluctuations of the field in terms of a 
measurement process [32], an aspect of this fascinating and successful theory which has been very 
little investigated until now. 
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